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Abstract 

The corrections to the current algebra mass formulae for the pseu- 
doscalar mesons are analyzed by means of a simultaneous exansion 
in powers of the light quark masses and powers of 1/N C . The rel- 
ative magnitude of the two expansion parameters is related to the 
mass ratio M^/M3, which represents a quantity of order N c m s . A 
set of mass formulae is derived, including an inequality, which leads 
to bounds for the ratios m u /m d and m s /m d . 
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The low energy properties of QCD are governed by an approximate, spon- 
taneously broken symmetry, which originates in the fact that three of the 
quarks happen to be light. If m u ,md,m s are turned off, the symmetry be- 
comes exact. The spectrum of the theory then contains eight strictly massless 
pseudoscalar mesons, the Goldstone bosons connected with the spontaneous 
symmetry breakdown. Their properties may be analyzed by means of an 
effective Lagrangian, which describes the low energy structure of the theory 
in terms of an expansion in powers of momenta and quark masses [0], 0. 

If the number of colours is taken large, the quark loop graph which gives 
rise to the anomaly in the divergence of the singlet axial current is suppressed. 
In the limit N c — > oo, QCD thus acquires an additional U(l) symmetry, whose 
spontaneous breakdown gives rise to a ninth Goldstone boson, the rf f|, [J. 
The implications for the effective Lagrangian are extensively discussed in 
the literature |3j and the leading terms in the expansion in powers of 1/N C 
have been worked out. The purpose of the present paper is to extend this 
analysis to first nonleading order and to discuss the consequences for the 
mass spectrum of the pseudoscalars. 

The relevant dynamical variable is the Goldstone field U(x), which lives 
on the quotient G/H, where G and H are the (approximate) symmetry groups 
of the Hamiltonian and of the ground state, respectively. Since I wish to 
study the large N c limit, G is the group 11(3)^x11(3)^ and £/(x)gU(3). The 
unimodular part of the field U(x) contains the degrees of freedom of the 
pseudoscalar octet, while the phase detU(x) = e 1 ^ ^ describes the rj'. 

The effective Lagrangian is formed with the field U(x) and its deriva- 
tives, C e ff = C e ff(U,dU,d 2 U, . . .). The low energy analysis relies on the 
symmetry properties of this function, which follow from the fact that the 
Lagrangian of QCD is invariant under independent U(3) rotations of the 
right- and lefthanded components of the quark fields - except for the U(l) 
anomaly which ruins the conservation of the axial singlet current and for the 
explicit symmetry breaking due to the quark masses. In an arbitrary chiral 
basis, the quark mass term is of the form q R mq L + q L m) ~q R . Using the pro- 
jectors P R = |(1 + 75), P L = |(1 — 75), the Dirac operator may be written 
as D = —ij^D^ + P L m + P R m) . The chiral rotation q R —>V R q R , q L ^ V L q L , 
V R , V L e U(3) takes this operator into D' = VD V, with V = P R V R + P L V L , 
V = 7°W7° = P L Vj + P R V^ . While the kinetic term is invariant under the 
operation, the mass term transforms according to m' = V R mV^. 

The path integral over the quark degrees of freedom is given by the deter- 
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minant of the Dirac operator. The singlet axial current is anomalous because 
the determinant does not remain the same under the above operation, but 
picks up a phase, det_D' = e iau detD, where a is the angle of the U(l) rotation, 
e ia = detV R V}, and 



In euclidean space, v is the winding number of the field configuration. Note 
that I am absorbing the coupling constant in the gluon field, = <9 M — iG^. 

The change in the Dirac determinant in effect adds a term to the La- 
grangian proportional to uj. So, if a term of this form is allowed for to start 
with, £qcd = £qcd — Ouj, the form of the Lagrangian remains the same under 
the full set of chiral transformations, but both the quark mass matrix and 
the vacuum angle 9 undergo a change, m^V R mV^, 9^ 9 — a. 

This property of QCD may readily be formulated at the level of the 
effective theory. The presence of an additional coupling constant, 9, also 
shows up in the effective Lagrangian, C e ff—>C e ff=£ e ff(U, dU, ... ,m,9). The 
expression must have the property that the transformation 



leaves it invariant, C e ff(U', dU', . . . , m', 9') =C e ff(U, dU, ... , m, #).[] 

Since the phase of the determinant = det U transforms according to 
0o = + a i the combination O + 9 remains invariant. It is convenient to 
replace the variable 9 by this combination, C e ff —£ e ff(U, dU, ... ,m, (fi + 9). 
The above symmetry relation then implies that the effective Lagrangian is 
invariant under a simultaneous rotation of the matrices U and m, at a fixed 
value of the last argument. The expansion in powers of derivatives and 
quark mass matrices yields a sequence of invariants which can be formed 
with m, U, dU, .... There are two differences compared to the corresponding 
series of terms which occur in the standard framework, where the rj' does 
not represent a dynamical variable: the determinant of U differs from one 
and the coefficients of the expansion are not constants, but functions of the 

More precisely, the corresponding action is invariant - the effective Lagrangian changes 
by a total derivative. The phenomenon gives rise to the Wess-Zumino-term. Since the 
corresponding vertices involve five or more meson fields, they do not contribute to the 
mass formulae and I therefore disregard the complication. 




u'=v R vvl 



w! = V R mVl , 9' = 9 -a 
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variable 0o + Up to and including terms with two derivatives or one factor 
of m, the most general expression consistent with the symmetries of the QCD 
Hamiltonian reads^ @, |§] 

Ceff = -V +V 1 (d fl tfd> 1 U}+V 2 (tfm}+V* (m^U)+V 3 d^d^ + ... , (1) 

with V n = V n (9 + 0o). The expression is valid to all orders in 1/N C . The 
symmetry does not constrain the form of the coefficient functions, except 
that Vo, Vi, V3 must be real and even under parity, V n (—x) —V n (x), while V 2 
obeys V 2 {— x) =^(2?)*. 

In the large N c limit, the coefficients reduce to constants. The reason is 
that, in this limit, the dependence of the various Green functions and matrix 
elements on 9 is suppressed 0: Compared to the leading contribution of the 
QCD Lagrangian, -G« v G a » v /Ag 2 , the term -9 G^G afiV /32n 2 represents a 
small perturbation of order g 2 6 ~ 0/N c . The vacuum energy density, e.g., 
is of the form e = N 2 g (9/N c ) + N c gi(9/N c ) + ... The ^-dependence of the 
coefficients occurring in the effective Lagrangian is similar 0. ||: 

K(0o + &) = (N C ) P "{V°(9) + N'W^e) + ...}, 6 = (0o + 0)/N c , 

with p = 2, pi =p 2 — 1,P3 — 0. In fact, the leading term in the expansion of Vo 
coincides with the gluonic contribution to the vacuum energy, Vq(x) = go(x). 

The above representation shows that the Taylor series of the coefficients 
with respect to (fio + 9 is accompanied by powers of 1/N C , so that only the 
first few terms are needed to a given order of the 1/N C expansion. The Taylor 
series of V , e.g., yields V ((f) + 6) = e + |r(0 o + Of + 0(1/N C ). The first 
term €q = 0(N 2 ) is an irrelevant cosmological constant, while r = 0(N®) 
represents the second derivative of the gluonic contribution to the vacuum 
energy with respect to 9 and is referred to as the topological susceptibility. 

In the present context, the dependence on 9 is of interest only to incorpo- 
rate the anomalous Ward identity and to control the expansion in powers of 
1/N C . The vacuum angle has now served its purpose and can be dismissed, 
setting = 0. Also, I now choose the chiral basis in such a manner that the 
quark mass matrix is diagonal and real. 

The preceding analysis involves a simultaneous expansion in powers of 
1/N C , powers of momenta p and powers of the quark mass matrix m. It is 

2 Note that the trace (d^UU^) coincides with id^o- 
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convenient to order this triple series by counting the three expansion param- 
eters as small quantities of order 1/N C = 0(5), p = 0(\J~5) and m = 0(5), 
respectively. Disregarding the cosmological constant mentioned above, the 
expansion then takes the form C e ff = C^J + C3 + . . . , where the first term 
is of order one,while the second collects the corrections of 0(5). Setting 
^i(O) = i-^ 2 ' ^2(0) = \F 2 B, the explicit expression for the leading term reads 

4$ = \F 2 (d^d»U) + \F 2 B(m(U + 17+)) - \r<g . 

The first two terms are familiar from the standard effective Lagrangian. They 
involve the pion decay constant F = 0(\/N2) and the constant B = 0(N®), 
which is related to the quark condensate. The third term is characteristic of 
the extension from SU(3) K xSU(3) Zy to U(3) fl xU(3) L . It equips the rj' with a 
mass proportional to the square root of the topological susceptibility. 

The following calculation also accounts for the corrections of first non- 
leading order. In addition to the terms arising from the Taylor coefficients 
Vj[(0) = 2iBK 1 and V 3 (0) = K 2 , which both are of 0(N?), these also in- 
volve contributions from higher orders of the derivative expansion, omitted 
in eq. ([!]). Their structure is known from the standard framework, where the 
corresponding effective coupling constants are denoted by L\, . . . ,L 8 . The 
first three of these do not play any role in the following, because they mul- 
tiply invariants of the type (dlf) A and do therefore not contribute to the 
masses of the pseudoscalars. The term Li(d^U^d^U)(mU^ + mil) contains 
two traces. Contributions of this structure can only arise from graphs with 
two or more quark loops: The term violates the Okubo-Iizuka-Zweig rule and 
is suppressed by one power of 1/N C . The same applies to the contributions 
proportional to L 6 and L 7 . The coupling constants L 5 and L s , however, are 
of 0(N C ) and do contribute to the masses at first nonleading order of the 
above expansion. Absorbing the constant B in the quark mass matrix with 
X = 2Bm, the effective Lagrangian becomes 

C eff = \F 2 {d^d»U + X (U + C/ f )) - Ho (2) 
+L 5 (dU^dU( X U + tf X )) + L 8 (U X U X + X^ X U j ) 
+iK 1 Mx(U' - U)) + K 2 d^ d^ + 0(5 2 ) . 

The first line contains the leading contributions of order N c p 2 ,N c m and 
7V C °, respectively. Their relative size depends on the relative magnitude of 
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the three expansion parameters 1/N c ,p,m. The second line contains the 
corrections of order N c p 2 m and N c m 2 , while the third one accounts for 
those of order iV c °m and N®p 2 . 

The mass spectrum is obtained by setting U = exp iip/F and working out 
the terms quadratic in the matrix field ip. In this notation, O is given by 
the trace (<p)/F, so that the quadratic terms are 

£eff = \{d^d^) + 2L 5 F- 2 ( X d^d^) + K 2 F~ 2 (d^)(d^) (3) 
- \{XV 2 ) - \tF- 2 (v) 2 - 2L 8 F~ 2 ( X pxV + *V> + 2K 1 F~ 2 (p)(xp) . 

For those fields which carry electric charge or strangeness, this expression is 
diagonal and yields 

M 2 + = (m u + m d )B{l + 8(m u + m d )(2L 8 - L 5 )BF~ 2 } , (4) 

and analogously for Mj^+ , M^ • These relations agree with those of chiral 
perturbation theory [0], with two simplifications: (i) The coupling constants 
L 4 and L 6 do not occur here, because they are suppressed by one power of 
N c . (ii) For the same reason, the chiral logarithms generated by the one loop 
graphs are absent. Since these graphs are inversely proportional to F 2 , they 
only show up at the next order of the expansion under consideration. 
The above mass formulae imply that the corrections in the two ratios^ 



Mjj m s + rh 



{1 + Am} 



Ml 2m 

m 

-{1 + Am} 



M^o - Mf + _m d ~m 
M 2 K - Ml ' m H -m 



are the same. Eliminating the quark masses in favour of M n , Mk, the explicit 
expression for the correction becomes 

A M = ^(M 2 K - M 2 )(2L 8 - L 5 ) . (5) 

In the double ratio 

Q 2^ M k M 2 k -Ml 



Ml M| - M 2 K+ 



The quantity rh denotes the mCctn mass of u cincl d 7 tH = \{jn u -j- tyi^ 



5 



the correction drops out. Accordingly, the corresponding ratio of quark 
masses is determined by the masses of the pseudoscalars 



^4^ = Q 2 » (6) 

m d~ m u 

up to and including first order corrections 0. Note that I have disregarded 
the electromagnetic interaction. Evaluating the corresponding self energies 
with the Dashen theorem [7[], one finds Q ~ 24. The theorem holds only in 
the chiral limit. The corrections of higher order are discussed in [||, |§. The 
corresponding uncertainty in the value of Q is of order 10 %. 

In the neutral sector, the contributions from L 5 and K 2 introduce off- 
diagonal elements into the kinetic term. These are removed with the change 
of variables (p—>-(p — 4{m, (p}L 5 BF~ 2 — 2((p)K 2 F~ 2 . The operation reduces 
the first line in eq.([5]) to the term ^(d^ipd^ip) and replaces the coefficients of 
the quadratic form in the second line by 

r = r(l-12K 2 F~ 2 ) , L 8 = L 8 - ±L 5 , K x = Ki + \K 2 + 2L 5 r F~ 2 . 

The masses of the neutral particles are thus obtained by diagonalizing this 
form. The mixing angles between the 7r° and T], rf are proportional to the 
isospin breaking mass difference m^ — m u . These angles are small, but play 
a crucial role e.g. for the transition 77 — > 3ir. In the masses of the neutral 
particles, however, isospin breaking only generates contributions of order 
(ma — m u ) 2 , which are negligibly small. Disregarding these, the ir° is degen- 
erate with tc^ and the 77 only mixes with the rj'. Setting (p = y^s^s + y?9V§; 
the quadratic part of the Lagrangian becomes 

£eff = K^sd'Vs + d^gd^ipg) - \{m\ipl - 20-1^9 + M 2 ipl) . 

The mass formulae for the charged particles may be used to eliminate the 
quark masses in favour of Mj^, M 2 . The coefficients then take the form 

m 2 = \{AM 2 K -Ml) + l(M 2 K -M 2 )A M 
a x = IV2(M 2 K -M 2 ){1 + A M - A N } 

M i= Q J~2 + |(2M|+M^)(1-2A JV ) + l(M 2 K -M 2 )A M . 
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Remarkably, m\ only involves the same correction Am which also occurs 
in the mass formulae for it and K. The term An = 12K\/F 2 , which de- 
scribes the Zweig rule violating contributions of order 1/N C , only affects the 
quantities o\ and M 2 . 

The eigenvalues M 2 = (M 2 , M 2 ,) obey (m 2 - M 2 )(M 2 - M 2 )=o\. Elim- 
inating Mi, this yields [m\ — M 2 )(M 2 , — m\) = o\ or, equivalently, 



Ml = < - ^ . (7) 

The relation states that the Gell-Mann-Okubo formula, M 2 = | (4Mj c — M 2 ), 
receives two categories of SU(3) breaking corrections: While the first is gov- 
erned by the same parameter Am which also determines the corrections in 
the masses of the charged particles and is accounted for in the term m 2 , the 
second arises from the rj-i]' transition matrix element of the operator qmq 
and is proportional to o\. 

The consequences of these mass formulae are the following. I first recall 
that the relation (||) constrains the two ratios x = m u /md and y = m s /md 
to an ellipse. Since Q 2 is very large, the ellipse is well approximated by 
x 2 + y 2 /Q 2 = 1, i.e. the center is at the origin and Q represents the large 
semiaxis, while the small one is equal to 1. To determine the individual 
ratios x, y one needs to know the quantity A M , which describes the strength 
of SU(3) breaking in the mass formulae and involves the two effective coupling 
constants L 5 ,L 8 . The former is known from the observed asymmetry in the 
decay constants F w , Fk, but, as pointed out by Kaplan and Manohar JTQ], L$ 



cannot be determined on purely phenomenologigal grounds. The scattering of 
the chiral perturbation theory results for the mass ratios m u /m^ and m s /m^ 
encountered in the literature [|ll|]-[[2(J] originates in this problem: Treating 
L 8 as a free parameter, one may obtain any value for Am and thus reach 
any point on the ellipse. In particular, the possibility that the mass of the 
lightest quark might vanish is widely discussed in the literature, because this 
would remove the strong CP problem p^ |. This possibility corresponds to 



rr iu/ m d = Qi fn s l v ^d = Q and requires that the correction is large and negative, 
A M = M 2 K /M 2 /(Q + \) - 1 ~ -0.45. 

I now wish to show that the framework specified above leads to a lower 
bound on Am which requires m u to be different from zero. The bound arises 
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because the relation (0) only admits a solution for <rn{, or 



AM 2 K ~ 3M„ 2 - Ml 

a « > - : { mi-md - -°-° 7 ■ < 8 » 

The inequality is an immediate consequence of the fact that 77-77' mixing leads 
to a repulsion of the two levels. Admittedly, the hypothesis that the first two 
terms of the 1/N C expansion yield a decent approximation for the theory of 
physical interest, N c = 3, goes beyond solid phenomenology This hypothesis, 
however, represents the only coherent explanation of the fact that the Okubo- 
Zweig-Iizuka-rule holds to a good approximation and I see no reason to doubt 
its reliability in the present case. 

The large negative value of Am required if m u were to vanish violates 
the above bound. The 1/N C expansion thus corroborates the conclusion that 
m u is different from zero. As is well-known, the quantity e %e det m then 
represents a physically significant parameter of the QCD Lagrangian. Why 
the mechanism which generates m and 9 happens to be such that, to a high 
degree of accuracy, this quantity is real, represents a well-known puzzle [^] . 
It cannot be solved the easy way, with m u = 0. 

A similar inequality also holds for the 1/N C correction, 

a ^ (y/3 + l)jg - (>/3 - 1)M* - 2MI 
N 8(M| - Mj) = °' 18 • 

This indicates that Zweig rule violating effects are not entirely negligible. 

As it stands, the mass formula ([^D only determines the magnitude of 
the SU(3) asymmetry Am~ (m s — rh) if the size of the Zweig rule violation 
An ~ 1/N C is known and vice versa. One may, e.g., choose a value for An, 
such that A M saturates the inequality This is the case, however, only 
if <Ji vanishes, i.e. if 1 + A N — A M = 0. In other words, the corrections 
would have to cancel the leading term. It is clear that in such a situation, 
the above formulae are meaningless. A coherent picture only results if both 
I Am I and \An\ are small compared to unity. The mass formula (0) shows 
that a negative value of the SU(3) asymmetry requires exorbitant 1/N C 
corrections. Even Am = calls for large Zweig rule violations, A^r — \- The 
term (1+A M — A N ) 2 , which corrects the contribution generated by the singlet- 
octet transitions for effects of first nonleading order, would then instead 
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modify this contribution by a factor of 4 if evaluated as it stands, eliminating 
it altogether if expanded to first order in the " corrections" . The condition 

A M >0 (9) 

thus represents a generous lower bound for the region where a truncated 
1/N C expansion leads to meaningful results. It states that the current algebra 
formula, which relates the quark mass ratio m s /m to the meson mass ratio 
Mk/M.1 represents an upper limit, m s /ih< 2M] C /M^ — 1. The corresponding 
bounds on the two ratios m u /md and m s /md depend on the value of Q, 

m u M*Q 2 -M*+ M\Ml 

m d M^ + M A K -M 2 K Ml ' 1 ] 

m s Ml(2Ml-Ml)Q* 

m d M$Q* + M A K -M 2 K Ml ' 1 ] 

If Q is taken from the Dashen theorem, these relations state that the Wein- 
berg ratios [12| correspond to the limiting case where the bounds are satu- 



rated, m u /m d > 0.55, m s /m<f< 20.1. Lowering the value of Q to Q = 22, the 
bound on m u /irtd becomes slightly weaker, m u /md> 0.48, while the one for 
m s /md decreases to m s /irid< 19.2. The estimate m u /m d = 0.3 ±0.1 obtained 
in ref . W? , |20|| from a multipole analysis of the transitions ip' — > ipir , ip' — > ipt] 



is outside this range and is thus not consistent with the above arguments (see 
also ref.0). 

To demonstrate that the observed mass pattern is perfectly consistent 
with the hypothesis that the corrections of order 1/N C as well as those of 
order m are small, I note that the contribution from the second term in 
eq.(0) is small, because it is suppressed by a factor of order MlfMl,. The 
corrections to that term are reduced by the same factor. The only new 
contribution in the second order formula which does not get suppressed is 
the one responsible for the difference between m\ and the value of which 
follows from the Gell-Mann-Okubo formula. Retaining only this term, the 
second order formula simplifies to 

4M| - Ml - 3M„ 2 2(M£ - M%) 



Am 4(M|-M2) + 3JWJ + Ml - AM 2 K ' 

Numerically, this gives Am = 0.18, thus requiring a breaking of SU(3) sym- 
metry in the mass formulae of the same order of magnitude as in the decay 
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constants. The corresponding value for the 1/N C correction is also quite small, 
An = 0.24. I repeat, however, that the above framework does not predict the 
values of the two corrections individually, but only correlates them. In par- 
ticular, the truncated large N c expansion is also consistent with a somewhat 
smaller breaking of SU(3) symmetry and a correspondingly larger violation 
of the Zweig rule. 
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